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7 = 2 Pion Scattering Length and Phase Shift with Wilson Fermions * 
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We present preliminary results of scattering length and phase shift for 1=2 S-wave tttv system with the Wil- 
son fermions in the quenched approximation. The finite size method presented by Liischer is employed, and 
calculations are carried out at /3 = 5.9 on a 24 3 x 60 and 32 3 x 60 lattice. 



1. Introduction 

Lattice calculations of scattering lengths and 
phase shifts for the pion system is an important 
step for understanding of strong interactions be- 
yond the hadron mass spectrum. 

There are already several studies of scatter- 
ing lengths in literature. For the 1 = 2 process 
calculations has been carried out with the stag- 
gered jl]|2| and the Wilson fcrmion actions 
f|. A recent work by JLQCD [Q carried out the 
calculation at several lattice cutoffs and obtained 
the value in the continuum limit. For the 1 = 
process, which is difficult due to the presence of 
disconnected contributions, a pioneering attempt 
has been made in Ref. Q . 

For the scattering phase shift, in contrast, there 
is only one calculation for / = 2 || , and the data 
are still too noisy to obtain reliable results. 

In this article we report on a high statistics 
calculation of the 1 = 2 S-wave pion scattering 
length and phase shift. We work in quenched lat- 
tice QCD employing the standard plaquette ac- 
tion for gluons and the Wilson fermion action for 
quarks. Simulations are made at (3 = 5.9 using 
200 and 173 configurations on a 24 3 x 60 and 
32 3 x 60 lattice. Quark propagators are solved 
with the Dirichlct boundary condition imposed in 



the time direction and the periodic boundary con- 
dition in the space directions. Quark masses are 
chosen to be the same as in the previous study of 
quenched hadron spectroscopy by CP-PACS || 
{i.e., m^/mp = 0.491,0.593,0.692,0.752 ) for 
both lattice sizes. 

2. Method 

The energy eigenvalue W p of an S-wave 7T7t sys- 
tem with momentum p and — p in a finite periodic 
box of a size 1? is shifted from twice the pion 
energy 2 ■ E p by finite-size effects. Liischer de- 
rived a relation between the energy shift AW P = 
W p — 2 ■ E p and the scattering phase shift S(p), 
which takes the form fft 



tan<5(p) 
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n ■ (27r/L) 2 , (n = 0,1,-- >„. 
AW p ■ (2E p L 2 )/(16tt 2 ) = 0(l/L), and v n , A n: 
and B n are geometrical constants. The scattering 
length is given by ciq = {ta,nS(p)/p) p ^Q. 

In order to obtain the energy eigenvalue W p 
we construct irir 4-point functions G p k(t) = 
(0\ilp(t)ttk(0)\Q) ■ Here £l p (t) is an interpolat- 
ing field for the S-wave 7T7t system at time t given 
by Q p (t) = Y, R n(R(p),t)ir(-R(p),t) where R is 
an element of the cubic group. In numerical cal- 
culations we construct the source operator f2fc(0) 
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Figure 1. R p (t) and D p (t) for p 2 = 0. 

using the noisy source method with U(l) random 
numbers. 

Since the 4-point function G p k{t) contains 
many exponential terms due to the Maiani- Testa 
no-go theorem @| , the extraction of energy eigen- 
values from G p k(t) is non-trivial. The following 
two methods are known |^| : 

(I) Diagonalization of G p k(t) at each t, where 
the momenta p and k are regarded as matrix 
indices. The eigenvalues take the form X p (t) = 
C p exp(-Wy) • [1 + 0(exp(-A p t)], where C p > 
and A p is the distance of W p from the other eigen- 
values Wk- In this method the fitting range of t 
should be taken large so that the 0(exp(— A p i)) 
terms can be ignored. 

(II) Diagonalization of M pk (t ,t) = [G{t )~ 1/2 
G(t)G(t )~ 1 / 2 ] p k at each t where to is fixed at 
some small value. The eigenvalues are given by 
X p (t) = exp(— W p (t - t )) where 0(exp(— A p t)) 
terms are absent. 

In the actual diagonalization we have to cut 
off the set of momenta. Here we expect that the 
components of G p k(t) or M p k{to,t) for p,k < q 
are dominant for the eigenvalue X q (t) in large t 
region, while the components p,k > q are less 
important. With this expectation, using both 
methods we calculate the energy eigenvalue for 
the ground state (n = 0) and the first excited 
state (n — 1) for V = 24 3 , and also that of the 
second excited state (n = 2) for V = 32 3 . The 
cut-off dependence is investigated by varying the 
number of momenta N > n for N = 0, 1, 2, 3. 
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Figure 2. R p (t) and D p (t) for p 2 = (2tt/L) 2 . 



3. Results 

In order to examine the effects of diagonaliza- 
tion, we calculate two ratios defined by R p (t) = 
G PP (t)/[G;(t)] 2 and D p (t) = A p (t)/[GJ (t)] 2 , 
where G p (t) is the pion propagator with momen- 
tum p. If G pp (t) or X p (t) behaves as a single expo- 



nential function, we can obtain the energy shift 
AW P easily from the ratio R p (t) or D p (t) by a 
single exponential fit. 

In Fig. [l] we compare the ratios for p 2 = 0. 
The two-pion source is placed at t = 8. The mo- 
mentum cut-off is set at N = 1 and both diag- 
onalizations are employed. The signals are very 
clear and diagonalizations do not affect the re- 
sult. We also checked the cut-off dependence by 
taking N = 2 and confirmed that it is negligible. 
In previous calculations of scattering lengths 
H the ratio Ro(t) was used to extract the energy 
shift A Wo- Our calculation demonstrates the re- 
liability of these calculations. 

We compare the ratios for p 2 = (2ir/L) 2 in 
Fig. ||. The momentum cutoff is set at N = 1 
and N = 2. The method (I) is used. In con- 
trast to the p 2 — case, the diagonalization is 
very effective. The cut-off dependence is negligi- 
ble, however. We also made the same analysis for 
method (II) and confirmed that the results are 
independent of the method. 

The analysis here, and additional one for p 2 = 
(2tt/L) 2 • 2 for V = 32 3 , lead us to conclude that 
(i) The momentum cut-off should be taken N > n 
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Figure 3. Results for scattering amplitude. 

for the energy shift AW P ( p 2 = (2tt/L) 2 -n ), and 
(ii) both method (I) and (II) give same results. 

In Fig. H we plot our results of the scattering 
amplitude T/(32tt) = E p ■ (tan£(p)/p) obtained 
by substituting our data for AW P into (0). The 
signals are clear in the small momentum region 
(upper figure in Fig. |||) but become noiser for 
larger momenta (lower figure). In order to obtain 
the scattering length and phase shift for various 
momenta at the physical pion mass, we extrapo- 
late our data with the following fitting assump- 
tion : T/(32tt) = A 10 ■ (ml) + A 20 ■ (ml) 2 + A 01 ■ 
(p 2 ) + A a2 • (p 2 ) 2 + An ■ (mlp 2 ). The fit curves 
are also plotted in Fig. [| 

In Fig. |^ we compare our results of the phase 
shift 8(p) at physical pion mass obtained with the 
fit above with experiments jl0],[ll]]. The simula- 
tion points are plotted by the large square sym- 
bols. Our results for 5(p) are 40% smaller in mag- 
nitude than those of experiments, and our result 
of scattering length, oom, = —0.0266(16), dif- 
fers from the ChPT prediction given by agm^ = 
-0.044 @. 

A possible cause of the discrepancy is finite 
lattice spacing effects. The JLQCD results for 
scattering length [Q show sizable scaling viola- 
tion, and hence that of the scattering phase shift 
cannot be considered small. Further calculations 
nearer to the continuum limit or calculations with 
improved actions are desirable. 
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Figure 4. Comparison of our results of scattering 
phase shift S(p) and experiments. 

10640248, 11640250, 11640294, 12014202, 
12304011, 12640253, 12740133, 13640260 ). VL 
is supported by the Research for Future Program 
of JSPS (No. JSPS-RFTF 97P01102). Simula- 
tions were performed on the parallel computer 
CP-PACS. 

REFERENCES 

1. S. Sharpe et al, Nucl. Phys. B383 (1992) 309. 

2. M. Fukugita et al, Phys. Rev. Lett. 71 (1993) 
2387; Phys. Rev. D52 (1995) 3003. 

3. R. Gupta et al, Phys. Rev. D48 (1993) 388. 

4. A. Aoki et.al, Nucl. Phys. (Proc. Suppl.) 83 
(2000) 83. 

5. H.R. Fiebig et.al, Nucl. Phys. (Proc. Suppl.) 
B73 (1999) 252 ; Few-Body System, 29 
(2000) 95. 

6. S. Aoki et.al, Phys. Rev. Lett. 84 (2000) 238. 

7. M. Liischer, Commun.Math.Phys. 104 (1986) 
177; 105 (1986) 135. 

8. L. Maiani and M. Testa, Phys. Lett. B245 
(1990) 585. 

9. M. Liischer and U. Wolf, Nucl. Phys. B339 
(1990) 222. 

10. W. Hoogland it et.al, Nucl. Phys. B 126 
(1977) 109 ; 

11. M.J. Losty it et.al, Nucl. Phys. B 69 (1974) 
185. 

12. J. Gasser and H. Leutwyler, Phys. Lett. B125 
(1983) 325. J. Bijinens et al, Phys. Lett, 
B374 (1996) 210. 



